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ABSTRACT 

Context. SAX J 1 8 1 8. 6-1703 is a flaring transient X-ray source serendipitously discovered by BeppoSAX in 1998 during an observa- 
tion of the Galactic centre. The source was identified as a High-Mass X-ray Binary with an OB SuperGiant companion. Displaying 
short and bright flares and an unusually very-low quiescent level implying intensity dynamical range as large as 10 3 ~ 4 , the source was 
classified as a Supergiant Fast X-ray Transient. 

Aims. The mechanism triggering the different temporal behaviour observed between the classical SGXBs and the recently discovered 
class of SFXTs is still debated. The discovery of long orbits (> 15 d) should help to discriminate between emission models and bring 
constraints. 

Methods. We analysed archival INTEGRAL data on SAX J1818. 6-1703. We built short- and long-term light curves and performed 
timing analysis in order to study the temporal behaviour of SAX J 1 8 1 8. 6-1703 on different time scales. 

Results. INTEGRAL revealed an unusually long orbital period of 30.0 ± 0.2 d and an elapsed accretion phase of ~ 6 d in the transient 
SGXB SAX J1818. 6-1703. This implies an elliptical orbit and constraints the possible supergiant spectral type between B0.5-1I with 
eccentricities e ~ 0.3 - 0.4 (for average fundamental parameters of supergiant stars). During the accretion phase, the source behaved 
like classical SGXBs. The huge variations of the observed X-ray flux can be explained through accretion of macro-clumps formed 
within the stellar wind. Our analysis strengthens the model which predicts that SFXTs behave as SGXBs but with different orbital 
parameters, thus different temporal behaviour. 
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1. Introduction 

Supergiant/X-ray binaries (SGXBs) are high-mass X-ray bina- 
ries (HMXB) composed of an OB supergiant companion star 
and a compact object (a neutron star (NS) or a black hole). In 
SGXBs, the compact object orbits within a few days (3-15 d) of 
the massive companion in circular (or slightly eccentric) orbits. 
Supergiant stars possess a strong stellar wind that is partly cap- 
tured by the compact object producing an X-ray radiation. This 
happens either by direct spherical accretion or by Roche-Lobe 
overflow via an accretion disk on the compact object. As the 
compact object always orbits within the stellar wind, these sys- 
tems are persistent X-ray emitting objects and show high varia- 
tions on short time scales. 

Very few of these systems were known before INTEGRAL'S 



launch ( Lutovinov et al.|2005 Dean et al.|20 05 ; Bodaghee et al.| 
2Q07\ . The mission allowed the discovery of several new sources 
that could b e identified as SGXBs (e.g. |Walter et al.|2006[|Chaty| 
et al. 2008). These sources show typical hard X-ray spectra of ac- 
creting pulsars and most of them show a strong absorption lead- 
ing to refer to them as obscured HMXB (e.g. [Filliatre & Chaty 
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* Discovery first reported in the 7th INTEGRAL workshop (Zurita 
Heras et al., in press). Based on observations with INTEGRAL, an 
ESA project with instruments and science data centre funded by ESA 
member states (especially the PI countries: Denmark, France, Germany, 
Italy, Switzerland, Spain), Czech Republic and Poland, and with the 
participation of Russia and the USA. 



2004; Rodrigu ezet al.|2006| |Bodaghee et al. 2006] |Zurita Heras 
et al. 120061 1 . Among the new SGXBs, INTEGRAL has unveiled 
a new subclass. Indeed, several newly discovered sources were 
identified as Galactic X-ray sources with a supergiant compan- 
ion and displaying a transient behaviour (e.g. IGR J 17544-2619, 
|in't Zand|2005[|Pellizza et al.|2006 1. The quiescent level in these 
systems is of the order 10^'" ergss " 1 (e.g. IGR J08408-4503 
|Gotz et al.|2007|[Leyder et al.|2007] l. The X-ray luminosity in- 
creases up to 10 JO ergss" 1 (as observed in known SGXBs) only 
during periods of short and luminous flares. The X-ray lumi- 
nosity remains at a very low level (if not totally absent) during 
most of the time, except during the flares. These flaring peri- 
ods last a few hours at most and then the source goes back to 
an undetectable level of emission (e.g. XTE J1739-302, Smith 



et al. 2006| l. Therefore, they received the name of Supergiant 
Fast X-ray Transient (SFXT, Negueruela et al. 2006[ ). Besides 
their transient nature, the spectral features of SFXTs are similar 
to the previously known persistent SGXBs. 

BeppoSAX discovered the new X-ray 
SAX J1818.6-1703 on March 11, 1998, during 



transient 
a bright 



flare ( |in 't Zand et al.| [l998]Q |in't Zand| < |2005) reported 



its possible association with other fast X-ray transients like 
XTE J1739-302 and IGR J17544-2619, archetype of SFXTs. 
Beside detections during deep surveys of the Galactic centre 
( |Revnivtsev et al.||2004] |Kuulkers et al.||2007) , similar short (a 
few hours) and bright flares of SAX J1818.6— 1703 were also 



All fluxes reported in the literature are summarised in Table [I] 
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Table 1. Previous detections summary of SAX J1818. 6-1703. 2. Observations and data analysis 



Energy range 


Flux 


phase 


keV 


mCrab 




2-9 1 


100 


0.24-0.28 


9-25 1 


400 


0.24-0.28 


18-60* 


4.9 


- 


18-45 2 


230 


0.10-0.13 


45-70 2 


170 


0.10-0.13 


20-30 3 


178 


0.12 


20-30 3 


185 


0.16 


0.5-10 4 


7.5 ±0.1* 


0.92-0.96 


20-60" 


11.1 ±0.9 




60-150" 


29.0 ±2.5 




18-45 5 


40 


0.06-0.08 


0.5-10 6 


< 0.11* 


0.56-0.59 


2-10 7 


4.2* 


0.13-0.17 



in unit of 10 12 erg cm 2 s_| 

BeppoSAX, 1998-03-11 ( |in 't Zand et al.|l998) 
INTEGRAL, Aug.-Sept. 2 003 JRevnivtsev et al. 1 2004) 
INTEGRAL, 2003-09-09 (Grebenev & Sunyaev 2005) 
INTEGRAL, 2003-10-09 and 10 (Sguera et al. 2 0lBj> 
Chandra, 2006-09-19, exp. 10 ks < |in't Zand et : al.|2006^ 
INTEGRAL, Feb.-A pr.+Aug-Oct' 2005 and Feb.-A.pr. 2006 
( |Kuulkers et al.|2007| ) 

INTEGRAL, 2008-04-15 and 16 dGrebenev & Sunyaev|2008) 
XMM-Newton, 2006-10-08, exp. 13 ks (Bozzo et al.|2008| ) 
Swift, 2008-04-18, exp. 2 ks jBozzo et aL)2008| > 



detected by INTEGRAL (see Table [T Grebenev & Sunyaev 



2005||Sguera et al |2005| [Grebenev & Sunyaev|2008| ). |Grebenev 



& Sunyaev ( |2005| l showed that the flares had a complex struc 



ture with a short 10-20 min precursor peak and a main 1.5-2 h 
long peak with different intensities. The hard X-ray spectrum 
(fitted with simple analytical models such as a power law) was 
also variable and became harder during the long flares. Using 
all INTEGRAL public data available before Mar. 2007, |Walter| 
|& ZuritaT -Ieras (2007 ) listed 11 flares with duration of 2-6 ks, 
and showed that SAX J1818. 6-1703 was always detected and 
displayed an intensity dynamic range of 250 between the lowest 
detection in deep mosaics and the brightest detection in indi- 
vidual pointings, in't Zand et al. (2006 ) and Bozzo et al. (2008 ) 



reported a quiescent low level of ~ (4 - 8) x 10~ erg cm" s 
with Chandra and Swift, and an unabsorbed 0.5-10 keV flux 
3cr upper limit of 1.1 x 10~ 13 ergcm^s -1 with XMM-Newton 
implying an intensity dynamic range of 10 3 ~ 4 as large as in other 
SFXTs. The Chandra spectrum was fitted with an absorbed 
power law: N H = (6.0 + 0.7) x 10 22 cirr 2 and F = 1.9 ± 0.3. 

Being similar to other SFXTs and searching a reddened 
early-type star in optical and infrared catalogues, Ne gueruela &| 
|Smith| ( |2006) proposed the source 2MASS J18183790-1702479 
as the likely counterpart. Confirmed by Chandra, the best 
X-ray position is R.A. (2000) = 18 h 18 m 37 9 s and Dec. = 
-17°0 2'47.9" (0.6" at 90% confidence level |in't Zand etaL] 
2006 ). The companion spectral type was constrained to a 09- 



Bl supergiant star (Neg ueruela & Schurch|200"7] ). 

In this study, we used the public INTEGRAL data in order to 
look for all the rare and short flares of SAX J1818.6-1703 and 
to study the long-term behaviour of the source, particularly its 
faint activity. In Sect. |2j we present the observations, the data 
analysis, the creation of deep mosaics and the timing analysis. 
In Sect. [3] we present the results and discuss the implications in 
SectH 



Launched in 2002 on a 3-days eccentric orbit, the ESA space 
mission INTErnational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL, Winkl eret al.|200 3 ) carries 4 instruments on-board 
of which 3 are dedicated to the detection of high-energy photons 
with good angular and spectroscopic resolutions and the last one 
is an optical V-band telescope. In this work, we used data from 
the hard X-ray and soft y-ray coded-mask imager IBIS/ISGRI 
whose main characteristics are a 15 keV-1 MeV energy range, a 
wide field of view (FOV) of 29° x 29°, and an angular resolution 
of 12' dUbertini et aLl200l||Lebrun et al.|2003) . 

The observations were divided in individual pointings with 
usual exposure time of ~ 1 - 5 ks. Discarding pointings with 
exposure times < 600 s and source off-axis angles > 14°, we 
analysed 4003 pointings distributed in 191 revolutions between 
Feb. 2003 (rev. (9046) and Nov. 2007 (rev. 0495). The total ex- 
posure time on the source was 9.3 Ms (i.e. 108.4 d) that were 
unequally distributed along the 4-years observations. 

The ISGRI data were reduced using the Offline Scientific 
Analysis version 7.0 software (OSA7). Individual sky images 
for each pointing were created in the 22-50 keV energy range. 
The source count rate was extracted in each image (pointings 
and mosaics) with the OSA tool mosaic_spec (version 1.7) 
where the position was fixed to the Chandra position and the 
point-spread function width was fixed to 6'. All detections of the 
source were considered above the 5cr significance level. Once 
the rare flares localised, we built deep-exposure sky mosaics 
where the pointings with detections were discarded. The goal 
was to observe the faint activity of the source and to constrain 
the lowest hard X-ray flux. We created mosaics for each revo- 
lution with exposures between 78-186 ks. Then, we combined 
consecutive revolutions in which the source was not detected 
and when the time gap between two revolutions was shorter than 
50 d. That resulted in 14 deep mosaics with exposure between 
200-1000 ks (except 2 cases with exposures of 60 and 80 ks). 
All the extracted fluxes are gathered to build the light curve 
of SAX J1818.6-1703 (see Fig. [T}. The 22-50 keV count rate 
can be converted with the relation 1 Crab = 177.5 counts s -1 = 



9.2x10 9 erg cm 
IJD + 51544. 



3. Results 



. The IJD time unit corresponds to MJD 



The source was detected in 32 individual pointings that were dis- 
tributed in 10 revolutions (crosses in Fig.[TJ). The count rate was 
mostly between (4.1 - 19.1) + 0.8 counts s -1 . On Sept. 9, 2003, 
there were 2 particular very bright flares that reach ed 44.8 ± 0.8 
and 39.1 + 1.0 counts s separated by 2 h (see 



Grebenev & 



Sunyaev 2005). Discarding all these pointings, we looked for 
the activity of the source around the flares within the same rev- 
olution. The source was only detected in 4 out of 10 revolutions 
with count rates between (0.9 ± 0.1) — (1.7 + 0.3) counts s _1 (tri- 
angles in Fig. [T] and 5cr upper limits when the source was not 
detected). In the other revolution mosaics, the source was only 
detected in 2 of them with similar count rates (boxes in Fig. [TJ. 
Both detections occurred before or after a revolution with flares. 
Beside these detections, we looked for other detections accu- 
mulating consecutive revolutions. Even with very long exposure 
time, the source was not detected with an average 5cr upper limit 
of 0.3 counts s" 1 (large 5o~ upper limits in Fig. [T}. In total, the 
source was detected 8 times for a few days at most. Detections 
1-2 and 3-4- were both separated by ~ 30 d, and 6-7 by ~ 60 d. 
The source seems to display a periodic recurrent flaring activity 



J.A. Zurita Heras and S. Chaty: The 30 days period of SAX J1818.6-1703 



3 



1 

+ 


i , , i , | 






1 




: t 
:.\ 

■ \ 

- t .h n 


< i.ni 



time [IJD] 

Fig. 1. 22-50 keV light curve of SAX J1818.6-1703. The data 
points are described in the text. 
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Fig. 2. 22-50 keV folded light curve of SAX J 18 18.6- 1703 with 
a period of 30.0+0.2 d and a time start of IJD 1 133.7365. Crosses 
correspond to ~ 2 ks flares, and boxes and 5cr upper-limits cor- 
respond to deep mosaics (typical exposure of ~ 100 ks). 



separated by ~ 30 d and lasting for ~ 5 d. However, flares in 
several observation periods were missed. 

To search for the missing flares, we studied the source 
long-term behaviour using the count rates extracted in point- 
ings and revolution mosaics. In pointings, the sensitivity dimin- 
ished when the source moved outside of the fully-coded FOV of 
ISGRI. The 5cr level varied between ~ 6 - 12 counts s _1 when 
the source off-axis angle varied between 9 - 14°. As the source 
was located at off-axis > 9° in ~ 2/3 of the pointings, the num- 
ber of flares below 12 counts s _I was underestimated. In mosaics 
with exposure < 20 ks, the 5<x level was usually > 4 counts s _1 , 
similar to the count rate of the faint detected flares. Therefore, 
the Lomb-Scargle periodogram (fast method of Press & Ry bicki| 



1989) was computed using the data of flares and mosaics with 
exposure > 20 ks. The best period was 30.0 ± 0.2 d (uncertainty 
derived using Eq. 14 of Home & Baliunas 1986 ). We folded the 
light curve with this period and a time start of IJD 1133.7365 
(see Fig. [2}. All detections fell between phases 0.0-0.2, so an 
elapsed time of 6 d. Only 31 out of 191 revolutions fell within 
the flaring periods. SAX J1818.6-1703 was detected in 12 of 
them as described above. We built a mosaic with the other 19 
revolutions (exp. 744 ks) and the source was detected with a 
count rate of 0.47 + 0.06 counts s , implying a variation fac- 
tor of 2-4 compared to the average count rates extracted in the 
6 revolutions with a detection. Combining revolutions per step 
of 0.2 phase between 0.2-1.0 (4 mosaics with exposure 1.5- 
2 Ms each), the source was not detected with a 5<x upper limit of 
~ 0.2 counts s _1 . 



4. Discussion 

As reported above, SAX J18 18.6-1703 was identified as a 
SGXB with a spectral type between 09I-B1I. Such systems 
usually radiate with an average luminosity of 10 36 ergs -1 with 
an important variability of factor < 20 and in which the com- 
pact object orbits around the supergiant companion at a distance 
< 2 R*. However, the source behaved differently compared to 
other known SGXBs. For SAX J1818.6-1703, we found a pe- 
riod of 30.0 + 0.2 d which can be interpreted as the orbital pe- 
riod of the system. A modulation of activity at a consistent pe- 
riod is found (D.M. Smith, private communication) in the RXTE 



Galactic bulge scan data ( |Markwardt et al. 2000| l. As observed 
by INTEGRAL, the source only radiated for ~ 6 d during its 
orbit. The detections reported by other missions are in good 
agreement with our ephemeris (see last column in Table TJ. 
During the 0.0-0.2 phase range, the average 22-50 keV flux 
in mosaics varied between (2 - 8) x 10 _1 1 erg crrT 2 s _1 . During 
short flares, the source usually reached fluxes of ~ (2 - 6) x 
10~ 10 ergcm^s -1 and, for very few flares, it even reached 
~ (1 - 2) x 10~ 9 ergcm^s -1 . During an XMM-Newton ob- 
servation corresponding to phase ~0.58 (i.e. near apastron pas- 
sage), the 0.5-10 keV 3<x upper-limit was < 10~ 13 ergcirT 2 s _1 . 
Therefore, the compact object must move on an eccentric orbit 
where it crosses the accretion zone only for 6 d. Only during this 
phase did SAX J1818.6-1703 behave as a SGXB. 

Let us consider the classical accretion model for SGXBs as- 
suming a smooth wind. The semi-major axis a of the system is 
determined with a 3 /P 2 = G{M„ + M x )/47r 2 with the period P, 
the stellar mass M„ and the mass of the compact object fixed 
to M x = 1.4 M Q . The stellar wind velocity of a supergiant can 
be approximated as v w j„d(r) = v 00 (l - R t /rf where is the 
terminal velocity, R* the stellar radius, and /3 = 0.8 - 1.2. The 
stellar mass-loss rate M wind = M wind (L t , M„ v^/v^, r eff ) is de- 



termined using eq. 12 of |Vink et al. (2000i where v esc is the 



stellar escape velocity, and r e ff the stellar effective temperature. 
Accretion can occur within a radius R ica - determined as R a 
2GM x /(v 



w nd + v x-* wnere v x is the compact object velocity. The 



r /4a 2 )M wind . 



accretion rate M accr is determined as M accr = (Rj 
The luminosity Lx can be approximated as Lx = eM accr c 2 with 
e ~ 0.2 the accretion efficiency. Considering typical stellar 
fundamental parameters of 09I-B1I stars (Martins et al.| 2005 



Crowfher et al. 2006} , the expected accretion rates vary only 



by factor 2-3 between spectral types. However, they strongly 
vary with the distance between the star and the compact object 
with factor ~ 10 between r = 2 - 3 R*. The expected fluxes 
when r — 3 R t are consistent with the low flux detected for 
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angle [degrees] 

Fig. 3. Top: SGXB with elliptical orbit (solid line), B0.5I star 
(solid circle), and e = 0.41. Circular dashed lines defined for r = 
2 and 3 R*. The accretion zone (r < 3 R*) is highlighted. Bottom: 
Solid line is the time of the compact object (left axis) along the 
orbit (f = d is aphelion); dot-dashed line is the distance star- 
compact object (right axis) along the orbit; and the light dashed 
lines determine the accretion zone where r < 3 R t and t e f accl . 



within the stellar wind. Indeed, it was shown that instabilities 
of the radiation line-driven outflows of hot stars generated struc- 
ture within the stellar wind affecting parameters such as the wind 
velocity, the outflow density, or the mass-loss rate (Runa cres &| 
Owocki 2002 and references therein). Recent studies of super- 
giant X-ray/UV lines emission have shown that optically-thick 
macro-clumps must form within the wind of hot supergiant stars 
( Oskinova et al.|2007| and references therein). Determining av- 
erage parameters of luminosity, frequency, and duration of flares 
in SFX Ts from INTEGRAL observations, |Walter & Zurita Heras| 
(2007) derived clump parameters that agreed with this macro- 
clump scenario developed by Oskinova et al. (2007), and pre- 
dicted that accretion in SGXBs should follow a two-regime be- 
haviour depending on whether the distance between the compact 
object and the supergiant companion was higher or lower than 
~ 2 - 3 R„ . |Negueruela et aT7| ( [2008| > also expected the same be- 
haviour and showed that the number of clumps that a compact 
object could statistically accrete within one orbit around a 08- 
B0I star varies between 1-8 with an orbital radius a — 3 - 6 R t . 
These values agree with the observation of SAX J 18 18. 6-1703 
flares, and can explain the observation of flares slightly outside 
the accretion phase as the flare observed by BeppoSAX. 

Acknowledgements. JAZH thanks Farid Rahoui, J. Rodriguez, D.M. Smith and 
R. Walter for useful discussions on SAX J18 18.6— 1703. The authors also ac- 
knowledge the use of NASA's Astrophysics Data System. 



SAX J 1 8 18. 6-1703 by the different missions. Therefore, we as- 
sumed that SAX J 1 8 18. 6-1703 starts to accrete when entering 
the zone within r < 3 R*. Still, the size of the accretion zone may 
vary up to 1 R* as the stellar mass-loss rate can change by a fac- 
tor of 2-3 in supergiants with identical spectral type ( |Crowther| 
|et al.|2006> . 

Under this condition and considering an elliptical orbij^) the 
period and the elapsed time of accretion f accr derived above im- 
plies constrains on the binary system eccentricity <?. For 091, 
09.51 and BOIa stars, the maximum elapsed accretion times al- 
lowed by the stellar parameters are 2.8, 3.0, and 4.6 d with 
e = 0.70, e = 0.69, and e = 0.58, respectively, that are in- 
ferior to the f accr derived above. For a B0.5Ia star, the con- 
ditions are fulfilled with e = 0.41 that implies a perihelion 
and aphelion of 2.3 and 5.5 R t (see Fig. [3j or e — 0.58 with 
a e 1.6 - 6.2 R*, respectively. For a Blla star, the conditions are 
fulfilled with e = 0.25 (i.e. a e 2.7 - 4.5 R t ) or e = 0.73 (i.e. 
a e 1.0-6.2 R t ). Only B0.5Ia and Blla stars seem to fill the 
accretion condition observed in SAX J 1 8 1 8. 6-1703. Discarding 
the short and bright flares detections, the X-ray flux varies be- 
tween ~ (2 - 8) x 10~" ergcm^s -1 along the whole ? accr . 
Therefore, the compact object distance r should not vary by a 
high factor when crossing the accretion zone, i.e. low eccentric- 
ities are favoured. Therefore, the case of a B0.5-lIa companion 
star with a low eccentricity (e = 0.41 or 0.25, respectively) are 
the best solutions to reconcile the wind model and observations. 
This determination is accurate to 0.5 spectral type since an ac- 
cretion zone of 4 R* would also give a solution for a B0I star 
with e = 0.16. 

Short and bright flares were also observed that were 10-100 
times higher than the average l evel during the ac creti on phase, 
and > 10 3 " 4 outside this phase. |in't Zand) | (|2005|> and | Walter & 



|Zurita Heras | ( |2007| l showed that such huge flares in transient 
SGXBs resulted from the accretion of massive clumps formed 



The Kepler's 2nd law allows to determine t = t(e, a, 9(r)). 
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